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Abstract 

We propose a Ginzburg-Landau model for the expansion of a dodecahedral 
viral capsid during infection or maturation. The capsid is described as a 
dodecahedron whose faces, meant to model rigid capsomers, are free to move 
independent of each other, and has therefore twelve degrees of freedom. We 
assume that the energy of the system is a function of the twelve variables with 
icosahedral symmetry. Using techniques of the theory of invariants, we ex- 
pand the energy as the sum of invariant polynomials up to fourth order, and 
classify its minima in dependence of the coefficients of the Ginzburg-Landau 
expansion. Possible conformational changes of the capsid correspond to sym- 
metry breaking of the equilibrium closed form. The results suggest that the 
only generic transition from the closed state leads to icosahedral expanded 
form. Our approach does not allow to study the expansion pathway, which 
is likely to be non-icosahedral. 

Keywords: 

2010 MSC: 92-XX, 20C40 



1. Introduction 

Most viruses are made of a protein shell, the capsid, built of identical 
protein units, that encapsidates and hence protects the nucleic acid (RNA 
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or DNA). Although viruses exhibit a wide diversity of shape a large number 
of them display icosahedral symmetry and this is irrespective of the number 
and the chemical composition of the protein subunits, the capsomers, that 
constitute the capsid. The basic principles to account for the icosahedral 
arrangement of the protein in a capsid were outlined by Caspar and Klug 
in the quasi-equivalence theory [Ti], and this continues to be a fundamental 
framework in virology. 

During their life-cycle viruses undergo structural transitions. These phe- 
nomena can be triggered by a change of the environment, such as a variation 
of temperature or pH. The occurence of a transitions induces a radial expan- 
sion of the capsid and a rearrangement of the capsomers and consequently 
the opening of pores on the capsid, so that the genetic material is exposed 
and eventually released in the host cell [2j [3J, IH [5] . 

In this work we focus on icosahedral viruses whose expansion can be mod- 
elled through the independent motion of twelve pentagonal blocks. These can 
be viral particles whose capsomers are pentamers (group of five proteins) or 
more complex viruses whose protein subunit are arranged to form pentagonal 
blocks [5]. 

In general, the energy of the capsid should account both for the cohesive 
forces between capsomers and for an internal pressure that tends to expand 
the capsid. However, instead of making specific choices, we adopt here the 
Ginzburg-Landau approach, and assume that the energy of the capsid is 
a function of the twelve independent degrees of freedom of the capsomers, 
and its explicit expression is given in terms of invariant polynomials of the 
icosahedral group. The coefficients of the expansion are the basic control pa- 
rameters of our model and, in turn, they should be related to the environment 
of the capsid. 

We relax the symmetry conditions and classify the minima of the energy 
according to their symmetry, in dependence of the parameters. The appear- 
ance of new minima corresponds to conformational changes of the capsid, 
with symmetry possibly lower than icosahedral. We restrict to minima of the 
energy corresponding to one of the three maximal subgroups of the icosahe- 
dral group: the tetrahedral group, and the dihedral groups of order 6 and 
10. 

The analysis shows that only four of the nine expansion parameters are 
relevant in this model: the parameter space is subdivided into regions, in 
which minima have a given symmetry. The main result of our work is that 
the only generic transition from a closed configuration with icosahedral sym- 
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metry turns out to be to an icosahedral expanded state. Minima with lower 
symmetry are not accessible from the closed reference state because they 
involve the cooperative change of more than one control parameters. 

Hence, our study shows that a phenomenological model based on the 
Ginzburg-Landau expansion is able to describe one of the main features of 
the expansion of a viral capsid: even though the transition pathway is unlikely 
to be icosahedral (cf. [El [7]), the final equilibrium state is still icosahedral 
(cf. 0). 

Our work complements other approaches to the study of conformational 
changes of viral capsids, either based on coarse-graining and phenomenolog- 
ical interaction potentials between the capsomers and relying on domain- 
decomposition techniques [HI [7] , or on the geometrical description of the fine 
features of the capsid via libraries of point sets with icosahedral symmetry 
[HI [TOl [HI [6] , normal mode analysis of the atomic ensemble of the capsid [12] 
or, finally, using the continuum theory of thin shells [131 EH]- 

2. Basics of the Ginzburg-Landau approach 

Let us consider a system described by state variables (xi, . . . , x n ) = x G 
W 1 and a symmetry group G acting on IR n . The action of G provides a 
representation p : G — > GL{R, n). We associate to the system a free energy 

E: R n x A — > E 

(x, cx) i — > E(x,a), 

where A C ]R m is an open subset, cx G A are parameters affecting the system 
(i.e., temperature, pH, etc.), and E G C 2 (lR n x A). We require the energy to 
be invariant with respect to p, i.e. 

E(p(g)x, a) = E(x, cx) Vg EG, Vx G R n Vc* G A. (2) 

The minima of E(-,a) with respect to x correspond to the stable phases of 
the system. Let us introduce the following definition 

Definition 2.1. Let p : G — > GL{R n ) be a representation of a finite group 
G. The isotropy subgroup of x G 1R™ is 

^ x = {g e G : p(g)x = x}. 
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Minima of E(-, a) can be classified according to their isotropy group. In 
particular, if Xq is a minimum of E(-,ct) such that 

we say that the system is in a high-symmetry phase. In general, we are 
interested in studying the local minima of the energy as ex. varies. The 
system undergoes a phase transition when the number and the symmetry 
of the minima changes as ex varies [To] . 

Notice that the level sets of E(-,a) are invariant under p. Therefore, if 
there exists a minimum Xq such that 

Saj = H, H < G, 

all points of the orbit p(G)xo also are minima. If G is finite, the orbit p(G)xq 
is also finite and has \G/H\ elements. Hence, low symmetry phases occur in 
different variants, while high symmetry phases are in general unique. 

Since we are interested in icosahedral viruses we consider the icosahedral 
group X which consists of all the rotations which leave a regular icosahedron 
invariant. It has order 60, and it is generated by the elements g 2 (twofold 
rotation) and g§ (fivefold rotation) such that (g 2 ) 2 = (fi's) 5 = (fi^s) 3 = e, the 
identity element. It is isomorphic to the alternating group of order 5, A$. 

Moreover X acts naturally on the twelve faces of a dodecahedron by per- 
mutation. In this way, we obtain a function o : X — > Si 2 , where Si 2 is the 
symmetric group of order 12, such that 

a(g 2 ) = (1, 6)(2, 5)(3, 9)(4, 10) (7, 12)(8, 11), 
a(^ 5 ) = (l,2,3,4,5)(7,8,9,10,ll). 

It is easy to see that a is an homomorphism, since 

°{92) 2 = (?{gsT = {a{g 2 )(x(g 5 )) 3 = id 5l2 . 

a is called a "permutation representation of X, being an homomorphism be- 
tween a group and a permutation group. The symmetric group S12 acts nat- 
urally on M 12 by permuting the indices of a vector. This action induces a rep- 
resentation p : X — > GL{R, 12) such that p(g 2 ) 2 = p(fi , s) 5 = (p(5'2)p(fi , 5)) 3 = 
-Z12, with J 12 the identity matrix of order 12. 
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3. Formulation of the model 

We now formulate a model for the expansion of the capsid in the frame- 
work of the Ginzburg-Landau theory. 

We focus here on viruses with capsid with icosahedral symmetry, whose 
expansion can be described via the independent motion of twelve pentagonal 
blocks. Therefore, before expansion, the capsid can be modelled eular 
dodecahedron (Figure [T]). We associate to each capsomer, labeled by i = 
1, ... ,12, a translation parameter X{ along the axis of the pentagonal face, 
with Xi G M (cf. Figure [2J. The translation variables x« have meaning only 
when non negative, but we allow them to assume negative values in order 
to simplify the analysis of the minima of the energy. However, we will only 
accept as physically meaningful the minima whose components are all non 
negative. 

The translation parameters Xi, % = 1, . . . , 12, represent the state variables 
of the system, i.e., x = (xi, . . . , Xi 2 ). Note that = (0, . . . , 0) represents the 




Figure 1: A dodecahedron: each pentagonal face is labelled by a number from 1 to 12. 
The face opposite to face i is labelled by i + 6, for i = 1, . . . , 6. 

closed configuration of the capsid, which has by definition full icosahedral 
symmetry, i.e., its isotropy subgroup is the whole group X. 

We require the energy E to be a C 2 function E : R 12 x A -> E, with 
AcR m , such that: 

(i) it is invariant with respect to the representation p of X on M 12 , i.e., 

E{p(g)x, a) = E(x, a) V</ G X, Va G A, \/x G M. 12 . 

(ii) = (0, . . . , 0) is a critical point of E, i.e. 

V x E(0, a) = Vck G A. 
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(iii) 

lim E(x, a) = +00 Va G A. 

\x\— >+oo 

As a first approach we also require the energy to be polynomial in the vari- 
ables Xi. This assumption relies on the fact that, if the energy is sufficently 
smooth, it can be expanded as a Taylor polynomial in a neighbourhood of 0, 
which still has the property of invariance with respect to p. 
Our goals are: 

• to find explicit expressions of the energy E as the sum of icosahedrally 
invariant polynomials; 

• to find the minima of E in dependence of the parameters a.. When 
0, which is always an extremum due to hypoyhesis (ii), is a minimum, 
the capsid is in the closed configuration. We interpret the appearance 
of new minima of the energy as conformational changes of the capsid 
associated to its expansion, driven by variations of the parameters ck, 
that in turn are related to the environment of the capsid. 




Figure 2: Basic variables of the model. 



4. Explicit form of the energy 

In this section we describe a method to find an explicit form of the energy. 
As we pointed out in the previous section, we assume it to be a polynomial 
invariant under the action of X. Therefore, we first study the ring of invariant 
polynomials and its structure. Complete proofs for the assertions below can 
be found, for example, in [16] or [T7] . 
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Let V be a vector space over a field K of characteristic zero (in what 
follows, we will be interested mainly in K = R) with dimV^ = n. We fix a 
basis A of V and let {x±, . . . , x n } be the dual basis of V*. 

Let G be a finite group and p : G — > GL(V) a representation of G. 
The action of G on V induces the dual action on V* and hence an action 
on the symmetric algebra Sym(V*) = K[xi, . . . , x n ]. We denote by K[x] G = 
K[xi, . . . ,x n ] G the subring of polynomials invariant under the action of G. 
K[x] is a graded ring, i.e., 

+00 

K[x] = ®K[x] k 

k=0 

where If [cc]fc is the vector space of the homogeneous polynomials of degree k. 
Since the action of G preserves the grading, K [x] G is a graded subring and 
we can write 

+00 

fc=0 

where K[x] G = K[xi, . . . ,x n ] G is the set of all the invariant polynomials of 
degree k. 

One way to find invariant polynomials is to take the average over the 
orbits of the action. More precisely, let us introduce the so-called Reynolds 
operator 

R G : K[x] — ► K [x] G 

defined by 

1 1 geG 

It is clear that Rq is a if-linear map that is the identity on the sub- 
ring K [x] G (in particular, it is surjective). Hilbert showed, as a consequence 
of his famous Basis Theorem, that 

Theorem 4.1 (Hilbert). The invariant ring K [x] G of a finite group G is 
finitely generated. 

Later E. Nother proved that the degree of the generators is bounded above 
by the order of the group: 
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Theorem 4.2 (E. Noether). Let \G\ be the order of the group G. Then 
K[x] G is generated, as a K-algebra, by (finitely many) polynomials of degree 
less or equal than \G\. 

Hence there exist p±, . . . ,pt £ K[x] G such that K[x] G = K\pi, . . . ,pt], 
i.e., every polynomial in K[x] G can be written as a polynomial in the pi 
and degpi < \G\ for i = 1, . . . , t. Moreover, we can find the generators by 
applying the Reynolds operator to a basis of the space of polynomials of 
degree less or equal than \G\, for example to the monomials. 

This gives the following 

Algorithm 4.1. In order to find a vector space basis of K[x] G , apply the 
Reynolds operator to all monomials in K[x]k. This yields a generating set of 
K[x] G as a vector space. By linear algebra, a basis can be extracted from it. 

To find a generating set for the full invariant ring, apply the previous step 
for k = 1,2,..., \G\. 

This algorithm has been performed successfully using the computer al- 
gebra software SINGULAR [18]. However, in the case where the order of G 
is large and the representation of G has a high degree, computations can be 
difficult, due to the high dimension of l(T[a;]&. One way to determine some in- 
variant polynomials with quite efficient computations is to use the irreducible 
representations of the group, as we are now going to discuss. 

Let p : G — > GL(V) be a representation of V , and let A be a basis of 
V , with dimU = n. Using Maschke's theorem, we know that there exist irre- 
ducible representations (defined over the algebraic closure of K) pi : G — > 
GL(Vi), i = l...r, with V = V x © . . . © V r , such that p = p x © . . . ffip r . Given 
a basis B{ of Vi, we have 



[p(g)]i 



( \pM\b^ ■■■ o \ 

V ... [p r (g)) Br J 



where B = B\ U . . . U B r and [p(<7)]s is the matrix of p(g) in the basis B. 
Moreover, if P is the matrix representing the change of basis from A to B, 
we have 

\p(9)\b = P- 1 \p(9)\aP V 9 eG. 

We denote by x' = (r)[ 1] , , r]i ] , VaLv^ ■••■> Vi\ ■ ■ ■ VdLv r ) and 
x = (xi, . . . ,x n ) the coordinates of a vector in W 1 in the bases B and A, 
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respectively, so that x = Px' . 

We consider the two rings of invariants K[x] G and K[x'] G . More precisely, 
p G K[x] G if 

p([p(g)] A x) =p(x) VgeG 

while / G K[x'f if 

f([p(g)] B x') = f(x') Vg G G. 

The two rings are clearly isomorphic, and an isomorphism is the following 

p:K[x'] G ^K[x] G 

f(x') f(x) = fiP-'x). 

In this way, it is possible to find polynomials of the ring K[x] G by first finding 
polynomials of the ring K[x'] G and then using the function tp. 

Let Bi = {v^}j, j = l...dimV^. We have, for every vector x' G lR n , 
since V = V 1 © ... © V r 

x' = x' 1 + ... + x' r 

with 

dimVi 
3=1 

Since [p] B = [pi]t3 1 © ... © [p r ]e r , invariance of / can be rewritten as 

/([pik^'i + . . . + \pAbX) = /K + ■■■ + <)• 

Therefore, we have 

($K[xr<ZK[x'\ G . 

i=i 

In other words, if pW(^) = p^( V ?, . . . , 77^) G K^, i.e. 

p^(\ Pl {g)] Bl <)=V {l \<) VgGG, 
then the polynomial 

£p ( V<) 



i=i 
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is an invariant polynomial in K[x'] G . This is a sufficient condition for the 
invariance which is useful for computations. 



Let us consider K[x'^, k > 0. Using Algorithm 4.1, it is possible to find 
asis {Phj}j of K[x'^\ G and write 



Since dimVi < dimV^ the computations are easier (sometimes in a consider- 
able way). We then consider the polynomial 

/(*')= E^E^te) 
k=i \i=i 

and finally using the function (p we find <p(f) = f(x) = f(P~ 1 x) e K[x] G . 
We organize the above results in the following algorithm 

Algorithm 4.2. Let V be a vector space with dimV = n, and A a basis 
of V. Let p : G — > GL(V) be a representation of a finite group G, and 
K[x\, . . . ,x n ] G = K[x] G the ring of the polynomials invariant under the ac- 
tion ofG. Perform the following steps: 

1. find the decomposition p\ © . . . © p r using the character table of the 
group G and the formula (see fiW ) 

Xp = X P 1 + ■ ■ ■ + Xprl ( 3 ) 

2. find explicitly the decomposition V = V\ © . . . © V r finding a basis Ei for 
each Vi using the projection operator (see JTffl ) Pi '■ V — > V,i defined by 

Pi = ^x Pi {g)p{g)i 

g&G 

3. find [pi\Bi, i = 1, ■ ■ ■ ,r and, letting B = B\ U . . . U B r , find the matrix 
P of the change of basis from A to B such that = P^pj^P; 

4. having fixed the degree d > and denoting by x' the coordinates of a 



vector in the basis B, find the polynomials pL-(a;Q forming a basis of 



K[x'j\^, using Algorithm 4-1, for i = 1, . . . r and k = 1, . . . , d; 
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5. write the polynomial 

k=\ \i=i \j=i / / 

where U^k = (\.im.K[x'^l and c*- G iT; 

6. perform the change of variables and find f(x) = f(P~ 1 x) E K[x] G . 

We now apply this algorithm to our model. In our case, G is the icosa- 
hedral group X, and p : X — > GL(R, 12) is the representation introduced in 
Section [2] Since X is isomorphic to the alternating group A 5 , its character 
table is the following 



Irrep 


C(e) 


C(g 5 ) 




C(92) 


C{g29b) 


Pi 


1 


1 


l 


1 


1 


P2 


3 


T 


l-r 


-1 





P3 


3 


1-r 


r 


-1 





Pi 


4 


-1 


-1 





1 


P5 


5 








1 


-1 



where r = + 2 , and C stands for conjugacy class. We note that, up to 
isomorphisms, there are 5 irreducible representations of X. We also note that 
since all characters are real, the direct sum decomposition is defined over R. 
In our case, we have Xpi.92) = and Xpi.95) — 2- We note that 1 + r + (1 — 
r) + = 2, and 1 - 1 - 1 + 1 = 0, so that 

Xp — Xpi + Xp2 ~l~ Xp3 Xp5' 

Using formula we find that p = pi © p 2 © ps © P5. Therefore, we have 
the decomposition 

R 12 = v 1 ®v 2 ®v 3 ® v 5 , 

where pi : X — > GL(Vi), and dimV2 = dimVs = 3, dimVi = 1, and dimVs = 
5. It is possible to compute the projection operators Pi, i = 1,2,3,5, and, 
as a consequence, a basis Bi = {v^ }j for each Vi, considering first the set 
{Pjej}^ 1; with A = {ej}j the canonical basis of IR 12 , and extracting a basis 
from it. We can then find [pijg, and the matrix P, whose columns are the 
vectors i = 1, 2, 3, 5, j = 1 . . .dimVJ, representing the change of basis 
from A to B = B x U B 2 U B 3 U B 5 . The results are 
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Representation p 1 



[Pl(52)]s! = [Pl(55)]s! = (1) 



• Representation p 2 



Ma*)] 



B 2 






-1 


-1 








-1 



Mas)] 



B-2 




Representation p 3 

[P3(#2)]s 3 = 




> [P3(5s)]s 3 




• Representation p 5 



Projection matrix 



/ 1 











\ 






f° 











-1 \ 














1 






1 











-1 











1 





1 


[ps(#5)k = 








1 





-1 








1 















1 








-1 




1 








) 






\o 








1 


- 1 / 



The following holds 



\M\b = P~ 1 \p{j9)\aP 



1 


-1 


Vs 


1 




-1 


1 


1 


1 


1 


-5 


1 


1 


-1 


1 


1 


-1 


Vs 


-1 


1 


1 


1 


1 


-5 


1 


-1 


-1 


-1 


1 


-1 




1 


1 


1 


1 


1 


1 


-1 


-1 


-Vs 


1 


1 


-1 


-5 


1 


1 


1 


1 


1 


-1 


1 


-1 


-1 


1 


1 


1 


-5 


1 


1 


1 


1 


-V5 


1 


-1 


-1 


-1 


-1 


1 


1 


-5 


1 


1 


1 


1 


-V5 


-1 


-Vb 


1 


-1 


1 


1 


1 


-5 


1 


1 


1 


-1 


-1 


1 


-VS 


1 


1 


1 


1 


1 


-5 


1 


1 


1 


1 


-1 


1 


-Vf, 


1 


1 


1 


1 


1 


1 


1 


1 


%/5 


-1 


-1 


1 


-5 


1 


1 


1 


1 


1 


1 


-1 


1 


1 


-1 


-1 


1 


-5 


1 


1 


1 


1 


V5 


-1 


1 


1 


1 


1 


1 


1 


-5 


1 


1 



/ [pi(#)k 



V 



B 2 



[pM] 



s 3 



\ 
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We are now able to find the invariant polynomials of degree k, pjy, form- 
ing a basis of M.[xi, . . . , Xdimvjf , using Algorithm 4.1 For the representation 



P2 and P3, which are 3-dimensional, we write the polynomials in the vari- 
ables x,y,z, and we perform the calculations for k = 2,4,6. There are no 
invariant polynomials of degree 1,3,5. For the representation ps, which is 
5-dimensional, we use the variables x, y, z, t, u. In this case we find invariant 
polynomials of degree 2,3,4. 



Representation p\ 

p ( H=x k k = l,2,... 

Representation p2 

k = 2 

(2)/ \ 2 2 2 2\/o , . 

P21 \ x i Vi z ) = x + y + z ^ — e - \ xz + v z ~ x v) 

k = 4 

(2)/ \ ( (2) x ~ 



Representation p 3 

k = 2 



jfc = 4 



pffV, V, z) = x 2 + y 2 + z 2 + (arz - yz - xy) 

5 



Representation p 5 

k = 2 



(5)/ , \ 2 1 2 1 2 1 j.2 1 2 

P21 \ x -iVi z,t,u) = x +y +z +t +u 

— (xy+xz + yz + xt + yt + zt + xu + yu + zu + tw) 
5 
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k = 3 

pfi (x, y, z, t, u) = xyz — xyt — xzt + yzt + xyu — xzu — yzu 
+ xtu — ytu + ztu 

pf^{x, y, z, t, u) = x 3 + y 3 + z 3 + t 3 + u 3 + -{xyt + xzt + xzu + yzu + ytu) 

5 

— -^( x2 y + xy 2 + x 2 z + y 2 z + xz 2 + yz 2 + x 2 t + y 2 t 

+ z 2 t + xt 2 + yt 2 + zt 2 + x 2 u + y 2 u 
+ z 2 u + t 2 u + xu 2 + yu 2 + zu 2 + tu 2 ) 

k = 4 

pfi(x, y, z, t, u) = x 2 y 2 + x 2 z 2 + y 2 z 2 + x 2 t 2 + z 2 t 2 + x 2 u 2 + t 2 u 2 
+ z 2 u 2 + y 2 u 2 + y 2 t 2 
+ xyzt + xyzu + xytu + xztu + yztu 

— - {x 2 yz + xy 2 z + xyz 2 + x 2 yt + xy 2 t + x 2 zt 

+ y 2 zt + xz 2 t + yz 2 t + xyt 2 + xzt 2 + yzt 2 + x 2 yu 
+ xy 2 u + x 2 zu + y 2 zu + xz 2 u + yz 2 u + x 2 tu + y 2 tu 
+ z 2 tu + xt 2 u + yt 2 u + zt 2 u + xyu 2 + xzu 2 + yzu 2 
+ xtu 2 + ytu 2 + ztu 2 ) 



p$(x, y, z, t, u) = x 4 + y 4 + z 4 + 1 4 + u 4 - -(x 3 y + xy 3 + x 3 z + y 3 z 

+ x 3 t + xz 3 + yz 3 + y 3 t + z 3 t + xt 3 + yt 3 + zt 3 
+ x 3 u + y 3 u + z 3 u + t 3 u + xu 3 + yu 3 + zu 3 + tu 3 ) 

+ -={x 2 yz + xy 2 z + xyz 2 + x 2 yt + xy 2 t + x 2 zt + y 2 zt 
5 

+ xz 2 t + yz 2 t + xyt 2 + xzt 2 + yzt 2 + x 2 yw + xy 2 u 
+ x 2 £-u + y 2 2-u + xz 2 u + yz 2 -u + x 2 tu + y 2 t-u + z 2 £u 

+ Xt 2 U + yt 2 M + Z 2 tu + Xt 2 M + yt 2 M + zt 2 U + XT/li 2 

+ xzu 2 + yz-u 2 + xtw 2 + ytu 2 + zt-u 2 ) 
[xyzt + xyz-u + xytu + xztw + yztiij 
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Let us now consider the basis B = {v^}^. Let us denote by 

a vector in IR 12 written in the basis B. Clearly, we have 
where 

dimVi 

i=i 

We can write an invariant polynomial F(x') G IR^'] 1 of degree d > as a 
sum 

fc=l \i=l,2,3,5 \ j / / 

The last step is to perform the change of coordinates from B to A. If we 
denote by ) the coordinates of a vector in IR 12 in the canonical 

basis A, we know that a; = Px'. Performing this change of variables in Q, 
we find 

fc=l \i=l,2,3,5 \ j 

where G IR. E^ d \x) is an invariant polynomial in IR[a;] :r of degree d > 0. 

We then choose d = 4 and, for semplicity, we discard the three poly- 
nomials of third degree p^ ,pi\ and p^ 2 ■ Further, in order for to be an 
extremum of the energy function, the linear terms p^(x) = Yl]=i x i mus t 
be absent. 

The output of the procedure is therefore a polynomial invariant energy 
of degree 4, which will be the focus of our subsequent analysis: 

E(x, a) =ap^ + bp^ + cp^i + dp?i 



(1) , (2) , (3) . (5) . (5) 



(5) 



where 

cx = (a, b, c, d, ci, C2, c 3 , c 4 , c 5 ) G IR 9 (6) 
and we assume that c, > 0, i = 1, . . . , 5 so that (iii) in Section [3] holds. 
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5. Bifurcation analysis 

In this section we characterize the minima of the energy in dependence of 
the parameters ex. The stationary points of the energy function are solutions 
of the system 

V x E{x, a) = 0. (7) 

We look for solutions of ^ with given nontrivial isotropy subgroup. The 
subgroups of I are listed in Table [TJ and the permutation representations of 
their generators are listed in Table [2] For a subgroup H of X, consider the 

Subgroup Generators Relations Order 



T 


92,g3d 


9l = 9 3 3d = (9293d? = e 
9ld = 9td = (95d92d) 2 = e 
9ld = 93 = {9392df = e 
9l = e 


12 


D 10 


92d, 95d 


10 


D 6 


92d, 93 


6 


Z 5 


95 


5 


D 4 


92d, 92 


gid = gl = {g2g2df = e 
gl = e 


4 


Z 3 


93 


3 


z 2 


92 


gl = e 


2 



Table 1: The subgroups of the icosahedral group. We have used the notation T for the 
tetrahedral group, Z?2n for the dihedral group of order 2n, and Z n for the cyclic group of 
order n. 

linear subspace of IR 12 whose elements are the vectors fixed by H: 

Fix(H) = {xe R 12 : p(g)x = x Vg £ H}. (8) 

Clearly, x £ Fix(H) if and only if if < Y, x . Therefore, in order to find 
minima of the energy with minimal symmetry H, it is enough to solve ((TJ) in 
Fix(if), i.e., letting 7 : R k — > Fix(if) C M 12 a parametrization of Fix(H), 
with k = dim(Fix(if)), to solve 

(V^o 7 )(t)=0, teR k . (9) 

This system has k independent equations only, since, omitting the depen- 
dence on cx for simplicity, the following holds 

(V^ o 7 )(t) e Fix(if), Vt G R k . 
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To see this, notice first that, since E is invariant with respect to p, then V X E 
is equivariant, i.e. 

V x E(p(g)x) = p{g)V x E{x) \/x G M 12 , 

so that, for t eR k , 

p(g)(V x E o 7 )(t) = p(g)V x E{<y{t)) = V x E(p(g)j(0) 
= V x E(j(t)) = (V x Eoj)(t), 

which is the assertion to be proved. 

Let a i = a(gi), . . . ,a r = a(g r ) be the permutation representations of 
the generators g\ . . . g r of H. Each <jj G Si 2 can be decomposed into dis- 
joint cycles, i.e. Oi = (nin 2 ■ ■ ■ n s )(mini2 • • • m r) ■ ■ ■ (P1P2 ■ ■ - Pt)- Let x = 
(xi, . . . , X12) G R 12 . Recalling that a { ■ x = (£^(1), . . . 2^(12)), if 

ry> ry> ry> 

•"fll • Lj U2 • • • ^TLs 

A mi • L m 2 ■ • • -^rrir . . 

< . (10) 

rf* rf rf 

x pi x p2 ... Xp, 

then <ji ■ x = x. Iterating this argument for all generators of the subgroup 
H we can find a parametrization of Fix(H). 



a(g 2 ) = (1,6)(2,5)(3,9)(4,10)(7,12)(8,11), 

a(g 2d ) = (1,12)(2,8)(3,4)(5,11)(6,7)(9,10), 

a(g 3 ) = (1,2,6)(3,5,10)(4,9,11)(7,8,12), 

a(g 3d ) = (1,10,2)(3,5,12)(4,8,7)(6,9,11), 

a(g 5 ) = (1,2,3,4,5)(7,8,9,10,11), 

a(g 5d ) = (1,10,11,3,6)(4,5,9,12,7). 



Table 2: Permutation representations and disjoint cycle decompositions of the generators 
of the subgroups of the icosahedral group. 

Table [3] lists Fix(if) for all subgroups of the icosahedral group. We only 
describe here the computations of Fix(Dio). In this case, the generators are 
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g 5 d and g 2 d, and (10) becomes 



x i — x w — x n — x z — x 6 

£4 = X5 = Xg = X\2 = x 7 



X X = Xn 

x 2 = X S 
X 3 = X A 

x 5 = xn 
Xg = Xio 



which combined yield 



X\ 
x 2 



x 3 
X 8 



X4 

y, 



x 5 — x 6 — x 7 — x 9 — x 10 — x ll — X 12 



which is the result in Table [3] 

For simplicity, we restrict our analysis to the maximal subgroups of X: 
the tetrahedral group T, and the dihedral groups D w and D e . 



Fix(T) 


= {x = 


(x, 


ry ry ry ry ry ry ry ry ry ry ry 
JU ■ JU ■ JU . JU - JU • JU ■ JU ■ JU - JU 2 JU ■ JU 


) : x e R}, 


Fix(L> 10 ) 


= {x = 


(x. 


1 1 If r f If If r f 1 1 If If If If 
y ■ JU ■ JU - JU - JU ■ JU y y ^ JU - JU - JU ■ JU 


:x,ye R}, 


Fix(L> 6 ) 


= {x = 


(x. 


x ,y,y,y, x , x , x ,y,y,y, x ) 


:x,ye R}, 


Fix(Z 5 ) 


= {x = 


(x. 


0C j tC j tC . 0C - %J - Z - Z . Z Z j Z j t J 


:x,y,z,te R}, 


Fix(£> 4 ) 


= {x = 


(x. 


y,z,z,y,x,x,y,z,z,y,x) 


:x,y,z<E R}, 


Fix(Z 3 ) 


= {x = 


(x. 


x ,y,z,y,x,t,t,z,y,z,t) : 


x,y,z,t G E}, 


Fix(Z 2 ) 


= {x = 


(x. 


y } z } t } y } x,u,w,z } t,w,u) 


: x, y, z, t,u,w G R}. 



Table 3: Fixed subspaces in M 12 of the subgroups of the icosahedral group. Notice that 
vectors in Fix(T) have full icosahedral symmetry. 



5.1. Stability of the closed capsid: conditions for the origin to be a minimum. 

Consider the energy ([5]). Since there are no terms of first degree, is 
always a stationary point of E for any a. 

The hessian matrix of E evaluated at has eigenvalues |a (multiplicity 1), 
Y^b (multiplicity 3), (multiplicity 3), and (multiplicity 5). Therefore, 

a, b, c, d > <^ is a minimum 

and the closed capsid is stable. 
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5.2. Solutions with icosahedral symmetry T or I 
In this case the system ^ reduces to 

1 1 „ 

-ax H — CiX = 0. 
6 3 

If a > it has only the trivial solution, while if a < (we always suppose 

Cj > for the positive defmiteness of the energy) we find x = We are 

only interested in strictly positive solutions, corresponding to the expansion 
of the capsid, so that we are left with the stationary point 



«. = «.-.«>. f = v / -^- (") 

The eigenvalues of the Hessian matrix H(xq) of E computed at Xq differ from 
those of H(0) only for the first one, which is now — |o. Therefore, when a < 0, 
loses its stability and a new minimum, Xq, appears, which corresponds to 
a complete expansion of the capsid which retains full icosahedral symmetry. 
The value of the energy computed at the icosahedral solution is 

E(x ,a) = -—. 

4ci 

5. 3. Solutions with fivefold symmetry D w 
In this case the system ^ reduces to 

^(5x + y) + A(y - x) + ^(5x + yf + ^(y - x) 3 = 0, 

i.e., 

(x - y) ( (x - yf H 

(5x + y) ^(5x + yf + ^r^j = 0. 

Assuming as usual that ci, c 5 > 0, for a > and d > the only solution is 
a; = y = 0. For a > and d < the solution 5x + y = is not acceptable 
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because one of the variables would be negative. For a < e d > the 



solution has icosahedral symmetry, and we recover the solution (11): 



y = x = 

For a < and d < the solutions are 



a 

2c x 



rp — I a |_ / d_ I ,,, / a / d_ 

X_,/ 2 C1 + V 2 C5 X - A/ 2ci V 2c 5 



(12) 



Both components of the first solution are positive for ajc\ < 25d/c^ < 0, and 
both components of the second solution are positive for ajc\ < djc^ < 0. 
Denoting by Xi,x 2 G M 12 the corresponding points, the eigenvalues of H{x\) 
are —j^d, — |a,^c (multiplicity 3), j^b ( multiplicity 3) and ^q(54 — 25C4/C5) 
(multiplicity 4). Therefore, if C4 > ffcs, ari< i c > 0, a?i is a minimum. 
Moreover, H{x2) has the same eigenvalues as H(x\), and therefore if C4 > 
|| C5 and 6, c > 0, also cc 2 is a minimum. 
To summarize, with the assumptions 

ad 54 
— < 25— < 0, &,c>0, ci,c 5 >0, c 4 > — c 5 , 
Ci c 5 25 

there are two stable phases X\ and x 2 after the expansion, whose isotropy 
subgroup is D 10 . The corresponding value of the energy is 



E(x x ) = E(x 2 ) 



c 5 a 2 + Cid 2 
4cic 5 



Instead, in the case where 



d a d 54 
25— < — < — <0, 6,c>0, ci,c 5 >0, c 4 > — c 5 , 
C5 Ci C5 25 

there is only one acceptable stable phase, corresponding to x 2 . 

Notice that for each solution considered here all other points in the icosa- 
hedral orbit are also solutions with the same symmetry. For instance, while 
x has icosahedral symmetry and therefore it is a full orbit, Xi and x 2 have 
symmetry D 10 , so that each corresponds to 60/10 — 1 =5 other equilibria 
with the same symmetry. 



20 




5.4- Solutions with threefold symmetry D 6 

We proceed as in the previous case. The system ^ reduces to 

x + y) + -^(x-y) + ^(x + y) 3 + -y) 3 = o 

x + y) (x — y) H -\x + y) (x — y) = 0, 

yj 60 v ' 24 v yj 1296 v ' 

i.e., 

(x + y) (^ + (x + y)'' 

, /108d , , 2 \ 

Assuming that c\, C4 > 0, for a > and d > the only solution is x — y = 0. 
For a > and d < the solution x + y = is not acceptable because one 
of the variables would be negative. For a < e d > the solution has 



icosahedral symmetry, as in (11), 

y = x = 

For a < and d < the solutions are 



a 

2c/ 



/ " 1 3 / 15d I a 3 / 15d 



2ci 1 5 V c 4 I ^ 4 A/ 2ci 5 V c 4 



_ / 3 / 15d I _ / a , 3 / 15d 

Both components of the first solution are positive for a/c% < ^-d/c^, and 
both components of the second solution are positive in the same range a/c\ < 
yd/ci. Denoting by x 3 , x± e IR 12 the corresponding points, the eigenvalues of 
H(x 3 ) are —j^d, —\a, (multiplicity 3), j^b (multiplicity 3), and j^d(25— 
54^) (multiplicity 4). Hence x 3 is a minimum if C5 > ||c4 and 6, c > 

Finally, the hessian matrix H(x±) has the same eigenvalues as H(x 3 ). 
Therefore, if C5 > ||c4 and b, c > 0, also #4 is a minimum. 

We conclude that, with the assumptions 

a ^ 54 d n i. n n 25 

— <- — < 0, 6,c>0, ci,c 4 >0, c 5 > — c 4 , 
Ci 5 C4 54 

there are two stable phases :r 3 and x^ after the expansion, whose isotropy 
subgroup is D§. The corresponding energy is 

E( X3 ) = E( X ,) = - C ^ph d L. 

4ciC 4 
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6. Conclusions 



We have explored in this work the possibility of using a Ginzburg-Landau 
approach to model conformational changes of viral capsids with icosahedral 
symmetry. Using a simplified model of the capsid as a dodecahedron whose 
faces, intended to model rigid capsomers, can move independently, we have 
assumed that the energy of the capsid is 4th degree polynomial invariant 
under the icosahedral group. The explicit form of the polynomial has been 
determined using the irreducible representations of the icosahedral group. 
The variables of the model are the amount of translation of each pentagonal 
face along its axis and, by consequence, we will only restrict to minima with 
nonnegative components. 

In this context, the basic parameters of the model are the coefficients 
of the polynomials corresponding to each irreducible representations, and 
the conformational change can be viewed as the change of structure and 
symmetry of the minima of the energy as the parameters vary. The phase 
diagram is depicted in Figure [3] when all parameters are positive, or when 
d < all other parameters being positive, the only minimum of the energy 
corresponds to the closed form of the capsid, but other minima appear, with 
non-icosahedral symmetry, in selected regions of the parameter space. 



a/c,= 25d/c 




Figure 3: Bifurcation diagrams for the minima of the energy for b, c, c\, C2, C3, C4, C5 > 0. 

It is believed that conformational changes in viral capsids occur via non- 
icosahedral pathways, but that the final form after maturation (of course 
when the capsid does not disassembly) is icosahedral. What is then the 
meaning of our non-icosahedral minima? We view the expansion process as 
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a conformational change of the closed, stable form of the capsid triggered by 
variations of the chemical environment of the virion. The destabilization of 
the capsid occurs by changing some or all the parameters in the energy func- 
tion, but, if the initial closed form corresponds to a point in parameter space 
in which all parameters are strictly positive, it is most unlikely that the min- 
ima cross special points of the boundary from the closed-form region to one 
of the low-symmetry regions. The most likely situation is that the minimum 
crosses a boundary, such as a = 0, and retains icosahedral symmetry. 

In other words, the conformational change leading to an icosahedral ex- 
panded state in viral capsids is preferred because it is generic, and therefore 
'easier': only one of the parametrs must change sign in order for an expanded 
icosahedrally symmetric state to appear. 
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